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MTIOML ADVISORY COMMITTEE FOB AERONAUTICS 
• AERODYNAMIC CHARACTERISTICS OF MCA RM-10 MISSILE 
IN 8- BY 6-FOOT SUPERSONIC WIND TUNNEL AT MACH 
NUMBERS FROM 1.49 TO 1.98 
III - ANALYSIS OF FORCE DISTRIBUTION AT ANGLE 
OF ATEACK (STABILIZING FINS REMOVED) 

By Roger W. Luidens and Paul C. Simon 


SUMMARY 

An analysis of the force distribution on a slender pointed 
body of rerolution at angle of attack -was made utilizing pressure- 
distribution data and balance measurements. The data were obtained 
in the MCA Lewis 8- by 6-foot supersonic tunnel at Mach numbers of 
1.49, 1.59, 1.78, and 1.98 and for a range of angles of attaok from 
0° to 9°. The Reynolds number based on the model length was approxi- 
mately 30,000,000. The parabolic body investigated was the half- 
scale model of the MCA supersonic flight-research' missile designated 
RM-10 (with stabilizing fins removed). A second model consisting of 
a cone-cylinder combination was investigated to isolate the effect of 
profile curvature. 

The inability to predict the normal force distribution due to 
angle of attack on slender bodies of revolution by the existing 
linearized potential theory was due in part to inaccurate prediction 
by the theory of the pressure distribution due to angle of attack on 
bodies with curved profiles, and in part to neglecting the effects 
of cross-flow separation by the theory. A concept of the linearized 
potential theory (in which the radius of the body is assumed to 
approach zero) is presented, which approximately eliminates the 
shortcomings of the theory with regard to the curved profiles. 

The axial friction and fore pressure force remained essentially 
constant with angle of attack. The increase in total axial force 
with angle of attack was primarily due to an increase in base pres- 
sure force. 
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INTRODUCTION 

Many methods are ‘ available for predicting the aerodynamic 
foroes acting on slender "bodies of revolution (at single of attack) 
moving at supersonic speeds. Two representative methods for pre- 
dicting the normal forces on slender todies of revolution inclined 
to a supersonic stream are linearized potential- theory, which seri- 
ously underestimates the measured values, and the method of ref- 
erence 1, which greatly improves the estimation ty modifying the 
potential theory in an attempt to account for viscous effects. The 
assumption of constant separation made in- reference 1, however, is 
inconsistent with .the pressure distributions observed on the RM-10 
model (reference 2). In addition, deviations of experimental pres-' 
sure distribution from potential theory exist over many regions of 
the body not appreciably influenced by viscous effects. 

An investigation was therefore conducted to determine the theo- 
retical and experimental distributions of the normal and axial forces 
and to study the origin of the discrepancies that occur between 
theory and experiment. Use was made of the RM-10 data presented In 
references 2 and 3 and the pressure-distribution data obtained with 
a second model. The investigation covers a range of Mach numbers 
from 1.49 to 2.00, angles of attaok from 0° to 9°, and a test 
Reynolds number of approximately 30,000,000 based on the length of 
the RM-10 model. 


APPARATUS AND PROCEDURE 

Data from the RM-10 (fig. 1(a)) investigation of references 2 
and 3 are presented herein. A second model (fig. 1(b)) having a 
length of 55 Inches, a maximum diameter of 2.5 inches, and a fine- 
ness ratio Z/D of 21.9 has also been investigated. The body con- 
sisted of a 12° vertex angle cone (nose was blunted by removing 
l/8 inch from the tip) extending 3.20 diameters, and a cylindrical 
section 15.4 body diameters long Joined by a curved section faired 
between the cone and the cylinder. Static-pressure orifices were 
-located In a longitudinal row on the model surface at every inch 
from station 2 inches to 17 Inohes and every 2 inches from sta- 
tion 17 inohes to 49 Inches. The support system described in ref- 
erence 2 was used for the cone-cylinder model. 

Experimental pressure data were obtained for nominal free- 
stream Ifeoh numbers of 1.5, 1.6, 1.8, and 2.0, for orifice radial 
positions of 9*0° and 180°, and for angles of attack of 0° and 
8° . The Reynolds number based on the model length was approxi- 
mately 23,000,000. 
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SYMBOLS 

The following symbols are used.- in this report: 
coefficient 


c a,c 


D 

f 

k 

M 

P 

q. 

S 

Sb 

S& 

u 


' x 
T0 


X,b 


x,r,0 


pressure coefficient, 


a constant 


P-PQ 

10 


section drag coefficient of a circular cylinder per unit 
length, in terms of its diameter 

maximum body diameter • 

local friction force 

a constant 

Mach, number 

static pressure - 

I ** 

Reynolds number, pllx/p. 
maximum cross-sectional area 
model base area 
model plan- form area 
Telocity 

radial yelocity component (cylindrical coordinates) 
axial velocity component 

tangential velocity component (cylindrical coordinates) 
coordinates of model 

cylindrical coordinates ( 9 =0 in plane of angle of attack 
and to windward) 
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a angle of attack 

$ cotangent of Mach. -angle, 'v/m 2 -! 

7 ratio of specific heats, 1.40 

A conditions for model due to angle of attack 

8 boundary- layer thickness 

5 * boundary- layer displacement thickness 

T) ratio of drag coefficient for circular cylinder of finite 

length to that for cylinder of infinite length 

H viscosity 

p density 

Subscripts : 

0 free- stream conditions 

2 conditions for model due to angle of attack 

A axial force (positive to rearward) 

a station axial force 

b base of model 

f due to friction 

L leeward side or surface of model 

2 length of model 

IT normal force (positive to leeward) 

n station normal force 

p pressure 

W windward side or surface of model 
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METHOD OF COMPUTATION AND SEDUCTION OF DATA 

The experimental pressure distribution and the force measure- 
ments for the EM- 10 model, which axe presented in references 2 and 
3, were reduced to normal and axial force coefficients. The incre- 
mental normal and axial pressure force coefficients (due to angle 
of attack) were obtained by graphically integrating the experimental 
incremental pressure coefficients C p 2 . Force coefficients are 

presented in terms of the model maximum cross-sectional area S with, 
the exception of the station force coefficients, which are based on 
the individual station diameters. 

Axial force data as a function of Mach number for angles of 
attack of 4° and 9° were obtained from cross plots because the 
experimental values recorded for the pressure model and the balance 
model were obtained at slightly different angles of attack. 

The theoretical incremental pressure coefficient associated 
with angle of attack for the linearized potential theory was given 
in references 2 and 3 as 


CL, o = 4a cos 0 — + a 2 (l - 4 sin 2 0) (l) 

dx 

Appropriate integrations (reference 3) of equations (1) over the 
body surface yield the following normal and axial force coefficients: 


40 n ,„ - *a g - 

a 2 

3 

(2) 

db 

AG n,L = rta a£ + 

a 2 

3 

(3) 

2 db 

AC a = -*a 2 5 


(4) 

s b 

Cn = 2a -g- 


(5) 

2 

AC a = -ar -g- 


(6) 
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The nor mal force coefficient was also computed by the formula given 
in reference 1: 


% ■ 



+ ^ c d. 



(7) 


Based on the conditions of this investigation and the criterions 
of reference 1, the value of tj was taken as 0.71 and c as 


1 . 2 . 


RESULTS AND DISCUSSION 

The aerodynamic forces acting on a body moving in a viscous 
fluid consist of the pressure forces associated with the acceler- 
ation of the fluid about the body and the shear forces resulting 
from the viscous attachment of the fluid to the body. For con- 
figurations with axial symmetry, these forces on the body are con- 
veniently analyzed in terms of components directed normal and par- 
allel to the body axis. 


Normal Forces 

A comparison of the measured total normal force coefficient 
and the values obtained by integration of the experimental pressure 
distribution coefficient with angle of attack for a range of Mach 
numbers is presented In figure 2. The difference between the total 
normal force and the normal pressure force Is the contribution of 
the friction, which is shown to be negligible. The conclusion that 
the normal pressure force represents the total normal force acting 
on the body is substantiated in part by the data- from two- . 
dimensional cylinders normal to a subsonic stream, which show that 
the measured total drag force agrees closely with the drag foroe 
calculated from the pressure distribution (for example, reference 4). 
Additional verification is presented in appendix A by an analytical 
analysis of the shear forces. The results of the analysis are 
included in figure 2 as zone A, which also evidences that the fric- 
tion mml force is very small within the angle of attaok range 
where the measured normal force is appreciable. (A small friction 
force does not preclude an effect of viscosity on the pressure dis- 
tribution; in fact this effe.ct may be appreciable ■, as indicated In 
reference 1.) 


* 
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Compared with, the experimentally observed variation of normal 
force coefficients Cjj with angle of attack a in figure 2 are 

the potential flow theory (equation (5)) and a modification of this 
theory given in appendix A, which accounts for the "boundary-layer 
displacement thickness observed at zero angle of attack. The 
increase in the normal force coefficient due to the boundary-layer 
displacement thickness is represented by zone B. Both theories 
yield a linear variation of Cjj with a‘ and do not approximate 

the experimental results, which lie considerably above these theo- 
ries and show an increasing normal force curve slope with increasing 
angle of attack. 

The method of reference 1 (equation (7)), which assumes a Vis- 
cous separation of the cross-flow along the length of the model, 
shows very close agreement with the data. The significance of this 
agreement, however, remains- to be deter mine d. None of the theories 
predict the increased normal force observed with increasing free- 
stream Maoh number.- . 

In order to facilitate further analysis of the normal pressure 
force distribution on the body, the experimental station normal foroes 
were separated into windward and leeward coefficients (AC n w and 

&c n,L> the s.um of which is the total station coefficient). The vari- 
ation of these coefficients with angle of attack for six represen- 
tative longitudinal model stations is presented in figure 3. For 
comparison, the local normal force coefficients calculated by poten- 
tial theory (equations (2) and (3)} are also presented. 

The potential theory best predicts the station normal force 
coefficient for both windward and leeward sides at the forward sta- 
tions. At the mid and aft stations the leeward normal force coef- 
ficients are greater than predicted by potential theory. This 
effect may not be arbitrarily charged to separation of the cross 
flow, as might be implied from the agreement shown in figure 2 
between the method of reference 1 and the total force. • The nor mal 
force coefficient, observed at the mid-section of the windward 
side of the body, which is not subject to separation, is also 
greater than that predicted by theory. The pronounced deviation 
of the exper ime ntal results from the theory on the aft-lee side of 
the body (stations 58 and 70) Is primarily associated with cross- 
flow separation. The agreement on the windward side in this region, 
however, is not adversely effected by the separation on the lee 
side. Thus, the possible argument that the cross-flow separation 
.of the leeward surface fore and mid-sections of the model may 
appreciably affect the normal force over the fore and mid-sections 
of the windward surface Is ruled out. 
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The excess lift, above that predicted by the theory, observed 
over both the lee and windward sides of the body mid-section results 
from an ant isymmetrical- type pressure distribution that contributes 
appreciably to the total body lift. The excess of lift observed 
ahead of station 45 inches is slightly less on the leeward than on 
the windward side, and it is noted that the effect of boundary- layer 
accretion (without the formation of vorticity) on the forward- lee 
side of the body (reference 2) is to decrease the lift on the lee 
side of the body. An insight into a possible origin of the pressure 
distribution that gives rise to the excess of lift observed over the 
mid-sections of the body may be found in the equation for the pres- 
sure distribution due to angle of attack (equation (1)), which con- 
tains a term that is antisymmetrical about the body (4a cos 0 

and results from the doublet distribution assumed on the body axis. 
The antisymmetrical term alone yields the normal force at a station. 
Thus it is logical to presume that the increased normal force, over 
that predicted by equation (l), observed on both the windward and 
leeward sides of the body mid-section arises for the most part from 
a shortcoming of the potential theory. The inability of the existing 
potential theory to predict the normal force coefficient therefore 
arises from (l) the failure of the existing theory to accurately 
predict the potential flow pressure distribution and (2) the neglect 
of viscous effects which result in cross-flow separation over some 
regions of the body. 

In an attempt to determine how the discrepancies between the 
experimental results and the potential theory which are not charge- 
able to viscous effects are affected by model contour, a pressure- 
distribution investigation was conducted on a second model (fig. 1(b)) 
composed of a cone gradually faired into a cylinder. The incremental 
longitudinal pressure coefficient distribution for the windward side, 
9=0° and the leeward side, 0 = 180° (the radial positions most 
significant with respeot to the pressure influence on the normal 
force) are presented in figure 4 for a = 8° at approximately the 
same Mach numbers as the EM- 10 data. 

The experimentally observed results show the best agreement 
with the potential theory on the windward side of the body for the 
regions where the body profile has constant slope (that is, the 
conical and cylindrical sections). On the lee side, close agreement 
was obtained over the cone but not over the cylinder, presumably 
because of cross-flow separation. Over the curved section of. the 
model, the experimental pressures at all Mach numbers sire greater 
on the windward and less on the leeward side of the model than those 
predicted by linearized theory. The same type of deviation from 
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theory was observed, in the EM- 10 data reproduced in figure 5. Pro- 
file curvature therefore has an effect on the pressure distribution 
that is not wholly accounted for by linearized potential theory. 
Appendix B presents a discussion of an assumption in the theory that 
results in this type of discrepancy and suggests a concept (which is 
not rigorous) that approximately corrects the shortc oming . 


Axial Force 

The axial force is the sum of the fore pressure force, the 
base pressure force, and the friction force. The distribution of 
the incremental station axial fore pressure force coefficient (due 
to angle of attack) AC & for several representative stations is 

shown in figure 6 as a function of angle of attack. Except for 
stations 6 and 18, the potential theory (equation (4)) accurately 
predicts the station axial fore pressure force on the body despite 
the presence of cross-flow separation and some shortcomings of the 
existing linearized potential theory, which have an adverse effect 
on the prediction of normal force. 

The increment of axial fore pressure force coefficient AC A p 

for the Mach numbers investigated is presented in figure 7 as a 
function of angle of attack. Also presented is the potential theory 
prediction of the axial force (equation (6)), which yields a thrust 
force proportional to the square of the angle of attack. The 
experimentally observed increment of axial fore pressure force coef- 
ficient is in general in the thrust direction, but is of much 
smaller magnitude than that predicted by theory. The high experi- 
mental axial fore pressure force as well as the experimental 
increase with increasing free-stream Mach number is associated with 
the flow about the forward portion of the body (fig. '6). It is 
emphasized that the maximum change in the Incremental axial fore 
pressure force coefficient shown in figure 7 represents a variation 
of 4 percent or less in the total axial force. This magnitude is 
evident in figure 8, which shows the total axial force coefficient 
C A and its component parts. Belative to the total axial force 

coefficient, the axial fore pressure force coefficient C A p 

decreases only slightly with angle of attack and is very nearly 
equal to the value predicted by linearized theory for a = 0 (ref- 
erence 3). Most of the increase in total’ axial force with angle of 
attack is accounted for by the increase in base force shown in 
figure 8(c). 
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The axial friotion force was calculated hy subtracting the 
axial fore and base pressure forces from the total axial force. 
Figure 8(d) shows that within the accuracy of the experiment the 
axial friction force coefficient C A f remains relatively con- 
stant with angle of attack, a result not readily anticipated in 
light of the pronounced circumferential variation in the boundary 
layer with angle of attack reported in reference 2. 


SUMMARY OF RESULTS 

An analysis of the force distribution over a slender pointed 
nose body of revolution at angles of attack between 0° and 9° was 
made from data obtained in an investigation conducted in the MCA 
Lewis 8- by 6-foot supersonic wind tunnel at free- stream Mach num- 
bers of 1.49, 1.59, 1.78, aid. 1.98 at a Reynolds number of approxi- 
mately 30,000,000. The following results were obtained: 

1. The contribution of the skin friotion to the total normal 
force observed on the body at angle of attack was negligible, and 
the normal pressure force represents the total force acting normal 
to the body axis. 

2. The existing linearized potential theory underestimates 

the normal force on the RM-10 for two reasons: (l) The theory does 

not accurately predict the potential flow pressure distribution due 
to angle of attack on bodies of revolution having curved profiles. 
(2) Separation of the cross flow exists over some regions of the 
body. 


3. The axial friction force and axial fore pressure force 
remained essentially .constant with angle of attaok. 

4. The increase in total axial force with angle of attack was 
primarily due to an increase in base pressure force. 


Lewis Flight Propulsion Laboratory, 

national Advisory Committee for Aeronautics, 
Cleveland, Ohio, 
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APPENDIX A 


BODY NORMAL FORCES AS AFFECTED BY 
FLUID VISCOSITY 

Viscous forces in a fluid may affect the normal force on a 
body of revolution inclined to a supersonic stream in several 
ways: (l) by shear forces acting on the body; (2) by changes in 

pressure on the body as a result of bo undar y- layer displacement 
thickness; and (3) by changes in pressure on the body due to flow 
separation. In this appendix the magnitude of the first two 
effects is examined. 

Fluid shear effect. - The contribution of the fluid Bhear 
force to the normal force depends on the local skin-friction coef- 
ficient and local inclination of the flow with respect to the body 
axis. The local friction force is 

df = CfCLQb d0 dx (8) 

The local inclination of the stream given in reference 2 is repre- 
sented by 


and 


Vg 

ss- = 2a sin 9 
u 0 


Jr 

*0 


db 

dx 


The station normal force may then be expressed as 

vrt pir 

rip = 2q. 0 b ( 2ct C f sin 2 9 d0 - C f cos 9 dSj 


O) 


( 10 ) 


Assuming that may be replaced by the constant mean value of 

0.0021 (based on wetted area) found experimentally to exist at 
Mq = 1.98 and a = 0, the station normal foroe coefficient per 
radian due to friction becomes 
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Integrating the station normal force over the length of the body 
yields, for the normal force coefficient due to friotion on the 
KM- 10, Cjj f/a «= 0.0764 (based on maximum cross-sectional area) 

compared with Cjj/a = 0.734 from potential theory, or approxi- 
mately 10 percent of the potential theory and 3 percent of the 
measured total normal force coefficient at a =* 8.61°, where the 
measured value is appreciable (fig. 2). It was found experimen- 
tally that the axial friction force remains constant with angle 
of attack, although the distribution of local friction undoubtedly 
changes. Such changes in local friotion about the body affeot the 
normal friction force through the first term in equation (10), but 
not appreciably. It may therefore be concluded that for the KM- 10 
the friotion normal force is small. 

Boundary- layer displacement thickness effect. - The equation 
for the normal force coefficient is 





cos 0b d0dx 


( 12 ) 


Utilizing linearized potential theory and including the boundary- 
layer displacement thioknees in the manner of reference 6, the 
equation for the pressure coefficient due to angle of attack 
(equation (l)) beoomes 

Cp g = 4a cos 0 + a, 2 (l - 4 sin 2 0) (13) 

This method assumes the boundary layer remains symmetrically dis- 
posed about the body at angle of attaok. For turbulent boundary 
layer of the form consistent with the profile observed on the RM-10, 
reference 2 gives for the boundary-layer thioknees 


and 


8 




(14) 


( 15 ) 
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Substituting in equation (13 ) and integrating in accordance with 
equation (12) over the EM- 10 where S = 1.0 inch, as observed in 

the plane of the base at M n = 1.98, k = 0.423 and ^ = 0.235 

u 5 

(reference 7), the normal force coefficient per radian due to the 
bou n dary- layer displacement thickness is §*/a = 0.216 as com- 
pared with the potential value of C^/a = 0.734. The contribution 

of the boundary layer to the normal force amo unts to approximately 
30 percent of the potential theory and 9 percent of the measured 
total normal force coefficient at a = 8.61°. 

The assumption that the boundary layer remains symmetrically 
disposed about the body at angle of attack is inconsistent with 
the experimentally observed results of reference 2, which showed 
a marked shift of the boundary layer from the windward to the lee- 
ward side of the 1 model. The present calculation, however, is 
useful to show an order of magnitude of the effect. 
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APPENDIX B 


EPFECT OP BODY PROPILE CURVATURE ON POTENTIAL 
PLOW PRESSURE DISTRIBUTION DUE TO 
ANGLE OP ATTACK 

Part of the discrepancy "between the experimentally observed 
and theoretically predicted pressures on the bodies investigated 
has been determined to be associated with the curvature of the 
profiles. It is now desirable to review the existing theory to 
determine what approximations in the theory may have led to this 
shortcoming and how the theory may be modified to account for this 
effeot. A comparison of the theoretical results, developed in ref- 
erences 2 and 5, for determining the flow about a body inclined to 
a supersonic stream reveals that the assumption of either subsonic 
or supersonic flow yields the same solutions. Beoause the develop- 
ment of the subsonic equations (reference 5) presents an easier 
physical picture of the quantities involved, this concept will be 
used. 


The equation for the perturbation velooity due to angle of 
attack is (reference 2) 

T x,2 * Q /<lb\ 

— 2a 008 e (s) 
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The implication in the subsonic case is that the component of flow 
normal to the body axis at any point on the body surface x^ of 
the preceding sketch is influenced by the body slope in the plane 
perpendicular to the body axis at x^. In the supersonio case, 

however , point (ii b-i 0 = 0) is influenced only by that portion 

j., j., 

of the body along and ahead of the fore Mach cone of the point 
under consideration. Because for the subsonic case the effect of 
angle of attack according to the theory arises only from the looal 
body slope in a plane normal to the axis through the point of 
interest, it is presumed that in the supersonic case the effect of 
angle of attaok arises only along the fore Mach line from the point 
under consideration. For bodies developed from profiles with cur- 
vature the fore Mach line traverses regions of varying slope 
(db/dx) between and Xg. The influence of the slope at vari- 

ous positions on the fore Mach line has not been determined; however, 
the use of an average value suoh as the slope at x^ - ^ pb-j_ was 

u 

used in this analysis. Other theoretical or physical averages could 
possibly be employed to advantage. It is of interest that suoh a 
concept does not influence the existing solutions in a cylinder or 
cone. The theoretical results yielded by this concept (which is 
not rigorous) are compared with the experimentally observed results 
in- figures 9 to 11 and show an improved agreement with experimental 
results except where cross-flow separation is known to exist. 


The distribution of incremental pressure coefficient as cal- 
culated by the stepwise doublet distribution method of reference 8, 


q 2 -2 7 ^»--, is presented in figure 9 for the 
u 0 


and defining Cp g 
EM-10 model when 0 = 0° and 180°, Mq = 1.98, and a = 8.61° 


This solution shows that in general the prediction of the magnitude 
and trends of the pressures is improved. The normal force coef- 
ficient obtained by the integration of this pressure distribution 
is presented in figure 11 and shows only a small improvement over 
the normal force determined by equation (5). 
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(b) Cone-cylinder model. 

' 

Figure 1. - Schematic diagram of NACA RJfl-10 half-scale model and cone 
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Figure 2* - Experimental and theoretical variation of normal force 
coefficient with angle of attaok at four Mach numbers. RM-10 
model. 
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(b) Leeward surface* 

Figure 3* - Experimental and theoretical variation of windward and leeward incremental station 
normal pressure force coefficient distribution with angle of attack for four Mach numbers* 
RM-10 model . 
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Figure 1. - Sjcperlmental and thaoretioal longitudinal variation of incremental pressure coefficient due to angle of 
attack for angle of attack of 8° and four Mach numbers. Cone-cylinder model; 0*0° and 180°. 
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(b) Leeward side of model; 6 * 180°.. 

Figure 5. - Bxperimental and theoretical longitudinal variation of incremental pressure 
coefficient due to angle of attack for -angle of attaok of 8.61° and three Mach numbers. 
Rif-10 model. 
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Figure 6. - Experimental and theoretical variation of incremental atation axial 
fore preaaure force coefficient distribution with angle of attack for four Mach 
numbers. RH-10 model. 
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Figure 7. - Experimental and theoretical variation of incremental axial fore pres- 
sure force coefficient with angle of attaok for four Mach numbers. RM-10 model. 




1423 









NACA KM E50I19 


23 


Ml 


Angle of 
attack, a 
(deg) 


(a) Total axial force coefficient. 


— Potential theory, a = 0 


(b) Axial fore pressure force coefficient. 



(c) Axial base pressure force coefficient. 
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(d) Axial skin friction force coefficient. 

Figure 8. - Variation of total and components of axial force coef- 
ficient with Mach number for angles of attack of 0°, 4°, and 9°. 
RM-10 model. 
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Figure 9 . - Experimental and theoretical longitudinal variation of incremental pressure 
coefficient due to angle of attack for angle of attack of 8.61° and three Mach numbers. 
RK-10 model. 
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Figure 10, - Experimental and theoretical longitudinal variation of incremental pressure coefficient due to angle of 
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Figure 11. - Comparison of modified theory with experimental normal 
force coefficient at angle of attack. RM-10 model. 
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